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ABSTRACT 

We use a statistically significant set of measurements to show that the field- aligned electron 
heat flux qu in the solar wind at 1 AU is consistent with the Spitzer-Harm collisional heat flux q s h 
for temperature gradient scales larger than a few mean free paths Lt > 3.5 A/ p . This represents 
about 65% of the measured data and corresponds primarily to high j3, weakly collisional plasma 
('slow solar wind'). In the more collisionless regime Xf p /Lr > 0.28, the electron heat flux is 
limited to q\\/qo ~ 0.3, independent of mean free path, where go is the 'free-streaming' value; the 
measured q|| does not achieve the full qo- This constraint q\\/qo ~ 0.3 might be attributed to wave- 
particle interactions, effects of an interplanetary electric potential, or inherent flux limitation. 
We also show a f3 e dependence to these results that is consistent with a local radial electron 
temperature profile T e ~ r~ a that is a function of the thermal electron beta a = a(0 e ) and 
that the j3 dependence of the collisionless regulation constraint is not obviously consistent with 
a whistler heat flux instability. It may be that the observed saturation of the measured heat flux 
is a simply a feature of collisional transport. We discuss the results in a broader astrophysical 
context. 

Subject headings: Sun: - solar wind, stars: winds, outflows 



1. Introduction 

Thermal conduction in the solar wind provides 
an important mode of energy transport and de- 
termines in part the radial electron temperature 
profile. The conductive, magnetic field-aligned 
electron heat flux is defined as 511 = — Ki|V||T e 
where K11 is the thermal conductivity coefficient 
and T e (r) is the electron temperature. In a fully 
collisional plasma, the Spitzer-Har m thermal con- 
ductivity ( Spitzer and Harm 19531 ) is km = ksh ~ 
3.16 n e T e T e /m e . Spitzer-Harm (SH) theory as- 
sumes that the electron distribution function f(v) 
remains approximately Maxwellian as it evolves 
through a temperature gradient scale, which cor- 
responds to assuming that the Knudsen num- 
ber K ~ Xfp/Lx is a small parameter , where 
L T = T e /\dT e /dr\ = R/a for T(r) - r" a , where 
R = 1 AU, and A/ p = v e r e is the mean free 
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path. In particular, SH theory assumes that the 
distribution function / rs fa + fi, with /o an 
isotropic Maxwellian and fi an anisotropic (cos 9) 
perturbation, such that ,fi < fp an d dropping 
0{fl) terms (jSpitzer and Harmlil953l) . The elec- 
tron collision time r e goes like Te^ 2 /n e , therefore 

5/2 

the thermal conductivity scales like ksh oc T e . 
Thus to maintain a constant conductive luminos- 
ity L — 47rr 2 <7|| , the wind must have T e (r) ~ r~ 2 / 7 . 
Measurements of the solar wind electron temper- 
ature profile generally show a power law profile 
Tpjr) ~ r~ a with valu es of a from 0.2 to 0.7 
(|e.g. Marsch et al.l 1989), not inconsistent with 
a = 2/7 rs 0.286. Of course, energy input from 
turbulent dissipation is likely to be important as 
well. 

The solar wind electron population at 1 AU 
consists primarily of a cool Maxwellian 'core' (~10 
eV, ~95% density), a suprathermal 'halo' (~70 
eV, ~4% density), and an antisunward 'strahl' 
population (-100-1000 eV, -1% density). The 
core is nearly isotropic, while the halo and strahl 
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exhibit clear temperature anisotropies. Since the 
electrons are subsonic (v e 3> v sw ), the maximum 
available heat flux corresponds to transport of the 
full thermal energy 3/2 n e kt,T e at the thermal 
speed v e ; this is often called the 'free-str eaming' 
or sa t uration heat flux q n = 3/2 n e ky,T e v e , (IParker 
1964t iRoxburghl Il97i ICowie and McKeefll977h . 
The ratio q s h/lo is proportional to the Knudsen 
number K ~ Xf p /Lr, the small parameter in 



SH t heory (jCowie and McKed Il977t I Salem et al 
20031) . 
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and this scaling of 311/50 with Xf p / Lt provides 
a clear test of SH theory. Assuming a value of 
a in Lt, we can compute the measured q\\/qo vs 
Xfp/Lr and compare to Equation (1) for q s h/qo] 
the only free parameter is a, the temperature gra- 
dient exponent. 

In this short paper we show that the measured 
q\\ in the solar wind is consistent with q s h until 
Xfp/Lx becomes as large at 0.3. Beyond that, in 
the collisionless regime, we find that 311 ~ 0.3 qo 
independent of mean free path. Collisionless reg- 
ulation of heat flux has been discussed in the 



conte xt of wave-particle interactions (|Garv et al 



1994 ) and escape from a n interplaneta ry electric 
potential (|Perkinsl [19731 iHollwed Il97i . We di- 
vide our data into intervals of electron thermal 
j3 e (ratio of electron thermal pressure to magnetic 
field pressure) and compare to theoretical thresh- 
old values for whistler and magnetosonic ins tabil- 
ity constraints on heat flux (|Garv et al 1 ll994l) . We 
find that the whistler instability overconstrains the 
measurements while the magnetosonic instability 
is more consistent. We also find that the data 
fit better to the SH relationship (in the collisional 
regime) if there is a /3 e dependence to the tem- 
perature profile index a. This may reflect the j3 e 
dependence of a wave heating mechanism, or may 
be a proxy for another parameter, such as colli- 
sional age. 

2. Measurements 

We use measurements of the solar wind elec- 
tron distribution function from th e Three Dimen- 
sional Plasma (3DP ) instrument ( Lin et al. 19951: 
Pulupa et all l2013t ) on the NASA Wind space- 



craft. The 3DP instrument uses two separate sen- 
sors - EESA-L and EESA-H - to measure the full 
3D distribution function f(v) in 88 angular bins 
from ~1 eV to ~30 keV, once per spacecraft spin 
(3 s). Eac h EESA sensor is a 't op hat' electrostatic 
analyzer ([Carlson et al.l Il982l ) designed to mea- 
sure solar wind thermal electrons (EESA-L) and 
suprathermals (EESA-H) each in 15 log-spaced en- 
ergy steps (few eV to 1.1 keV for EESA-L and 100 
eV to 30 keV for EESA-H). Measurements from 
the two detectors are combined to form a single 
distribution function and this distribution func- 
tion is corrected for spacecraft floating potential 
effects using quasi-thermal noise measurements as 
an absolute plasma density benchmark; low en- 
ergy monopole corrections (few Volts) are less im- 
portant for higher-order moments such as heat 
flux, however dipole fields ma y introduce errors o f 
5% to the heat flux moment ( Pulupa et al.ll2013h . 
We use 155,182 independent measurements of f(v) 
from two, 2-year intervals: 1995-1997 (solar min- 
imum) and 2001-2002 (solar maximum) and in- 
clude only 'ambient' solar wind (no CMES, fore- 
shock, etc.). Intervals of 'bi-directional' heat flux 
(usually associated with CMEs) are also excluded. 
We compute the electron heat flux moment as 



q e = - j '/>"/( c'\ir-ir 



(2) 



from the measurements, where the secular ve- 
locity is w — v — Vf, and Vb is the bulk speed. Here 
we consider the (dominant) magnetic field-aligned 
component of the heat flux gn = q-B. The thermal 
electron properties qo, Xf p , and j5 e are computed 
using measurements of the 'core' electron popula- 
tion, a point that we discuss below. 

Figure 1 shows the joint probability distribu- 
tion of q\\/qo and Xf p /Lx normalized to the peak 
value in each Xf p / Lt histogram and Equation (1) 
is over plotted as a diagonal line in the top panel. 
The number of points in each A f p /Lx bin is shown 
in the lower panel and a temperature exponent of 
a ~ 2/7 is used to calculate Lt- It is apparent 
that q\\/qo tracks Equation (1) over much of its 
range. 

Figure 2 shows the mode (most probable value) 
in each bin of A/ p /Lt- The modes are calculated 
directly from the distribution of data (red points) 
and from log- normal fits (black points). In this 
Figure and in Figure 1, the Spitzer-Harm rela- 
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Fig. 1. — Joint normalized distribution of the nor- 
malized electron heat flux q\\/qo and the temper- 
ature Knudsen number K — Xf p /Lr (top panel) 
and the distribution of data with X/ p /Lt (bottom 
panel). The diagonal line in the top panel is the 
Spitzer-Harm relationship (Equation (1)). 



tionship in Equation (1) appears to be a reason- 
ably good approximation to the data. At about 
Xf p ~ 0.28 Lt, the measured data breaks from the 
SH line and flattens to a fixed value of g|| ~ 0.29 go- 
Approximately 65% of our measurements lie in 
the SH regime, the other 35% in the collisionless 
regime g|| ~ 0.29 go (bottom panel of Figure 2). 
The so-called 'collisional age' A = v . w t is often 
used to measure collisi onal evolution ( Salem et al.l 
20031 iBale et al.ll2009h . especially when consider- 
ing interactions with ions convecting at the so- 
lar wind speed v sw . The collisional age correlates 
strongly with solar wind speed; fast solar wind is 
more collisionless (hot and rarified) and slow wind 
is more collisional (cooler and dense); therefore 
data in Figure 1 and 2 in the collisional regime 
(Xfp/Lx <C 1) is primarily slow wind, while colli- 
sionless data (Xjp/Lt 1) is primarily fast wind. 

3. Electron (3 dependence 

In Figure 3, we break this data into 4 intervals 
of electron thermal beta j3 e = (n e kBT e )/ (B 2 /2/iq). 
Since Xf p /Lx oc l/n e and density variations dom- 
inate the pressure variations, the high /3 e plasma 
corresponds to collisional plasma (small Xf p / Lt)- 
This evolution can be seen in the panels of Fig- 
ure 3, organized from high to low /3 e (top to bot- 
tom); the points move towards the right (towards 
the collisionless regime). As the clusters of points 
move to the right, they maintain the SH-like power 
law behavior, but require different values of a (the 
temperature profile index) to conform to the curve 
- note that a is the only free parameter. Alignment 
to the SH curve give corrected values and suggest a 
j3 e dependence to the electron temperature profile 
index a and to the breakpoint A& between the col- 
lisional SH and the collisionless heat flux regimes. 
Table 1 summarizes these results. 

Figure 3 also shows the /3 e -dependent instability- 
limited electron heat flux values for both the mag- 
netosonic and whistler instabilities as calculated 
by Gary et al. (1994). We use most probable 
values of j3 e in each interval and scale the corre- 
sponding (jmax = 10~ 3 Q p ) threshold from Figure 
1 of Gary et al. (1994): magnetosonic-limited val- 
ues are shown as horizontal dotted red lines, while 
whistler instability-limited heat flux is shown as 
blue dotted lines. It can be seen that the heat flux- 
driven whistler instability overconstrains the data; 
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Fig. 2. — Plot of the most probable values (modes) 
of q\\/qa in bins of Xf p /Lx with standard error er- 
rorbars (very small) . Red points are the computed 
modes, black points are modes determined from 
log normal fits. Again, the black diagonal line 
in the upper panel is the SH Equation (1). The 
measured values depart from SH behavior above 
Xfp/Lx ~ <?||/<7o ^ 0.28. The cumulative distribu- 
tion in the bottom panel shows that approximately 
65% of this dataset corresponds to SH heat flux. A 
symbol at qu — qo show the free-streaming value. 




Fig. 3. — Data from Figure 2 broken into four in- 
tervals of electron thermal /3 e . In each interval, the 
data are consistent with the SH relationship with 
the temperature gradient index a as a free param- 
eter. The /3 e dependence of a is summarized in 
Table 1 . Dotted horizontal lines are the thresholds 
of the magnetosonic instability (red) and whistler 
instability (blue), both from Gary et al (1994). 
The whistler instability appears inconsistent with 
the heat flux levels in the collisionless regime. 
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the magnetosonic instability underconstrains the 
data in the SH regime, but is in fact consistent 
with the limiting electron heat flux in the colli- 
sionless regime. 





N 


a 


A 6 


0.01 - 0.4 


12,507 


0.51 


0.40 


0.4 - 1.5 


50,056 


0.31 


0.27 


1.5 - 5.0 


80,560 


0.20 


0.18 


5 - 100 


12,059 


0.13 





Table 1: The (3 e dependence of T e (r) ~ r~ a and 
the SH-collisionless transition scale K^Lt- N is 
the number of measurements in each f3 e inter- 
val. Note that the limiting heat flux in the col- 
lisionless regime is proportional to the breakpoint 
q\\ ~ 1.07 A fc g - Also note that 2/7 « 0.286. 



4. Conclusions 

In this manuscript, we show that the transition 
from collisional (Spitzer-Harm) to collisionless 
electron heat conduction in the solar wind occurs 
at a mean free path of about A/ p ~ 0.28 Lt, where 
Lt is the electron temperature gradient scale. In 
the collisionless regime, the heat flux is limited to 
<7|| ~ 0.3 go, where go is the free-streaming value. 
Some previous analyses suggested departures from 
the S H value (jFeldman et al.l Il975c iPilipp et al 



19871 ). This could be attributed to the choice of 
data intervals in those analyses, none of which 
were very statistical. In fact, a preliminary statis- 
tical analysis of the Helios electro heat flux mea- 
surements shows similar results to what we present 
here (K. Horaites, private communication). 

In our analysis, we have used the 'core' electron 
density and temperature to compute go and Xf p . 
Theoretical work on heat conduction on steep tem- 
perature gradients has shown that a self-consistent 
flux limitation arises as SH theory begins to break 
down (for /i ~ / ) a nd that this correspo nds to 



order Xj p /Lt ~ 0.1 (IShvarts et al.l 119811 ). simi- 
lar to our results. The growing departure of the 
measurements from the SH curve in our Figure 2 
is consistent with this effect (which can be seen 
by multiplying Figure 2 of (jShvarts et al. 1981 ) 
by Xfp/Lr)- Since / ~ n/T 3 / 2 , the constraint 
/i < fo is well-satisfied for solar wind parame- 
ters (n h ~ 1/20 n c and T h ~ 7 - 10 T c ) for ve- 



locities less than about v < 2.6 v t h. c , therefore 
the core population does represent the collisional 
physics. Similar results were obtained by Smith et 
al. (2012), who solved the electron kinetic equa- 
tion with a linearized Fokker-Planck operator in a 
fixed ion profile and found limiting heat flux values 
comparable to ours and attribute it to skewness in 
the distribution function at speeds of v > 3 v t h- A 
simple model of a suprathcrmal tail escaping from 
an interplanetary electric field (lHollwedll97"4l) pre- 



dicts a collisionless heatflux <; ~ 3/2 n e k^,T e v sw , 
which is also approximately consistent with our 
results, since v sw ~ 1/3 v e . If we repeat our 
analysis with the full electron temperature T a u = 
(n c T c + n h T h + n s T s )/(n c + n h + n s ), we find a 
similar breakpoint, but somewhat faster departure 
from SH. 

Our measurements show that the solar wind 
heat flux is well-described by the collisional 
Spitzer-Harm value until the mean free path is ap- 
proximately one third of the temperature gradient 
scale Xf p ~ 1/3 Lt and for larger A/ p is propor- 
tional to the saturation value q ~ 1/3 go- While 
the limiting mechanism is not yet understood, 
these results should be useful for solar wind and 
coronal modeling efforts. Our analysis also sug- 
gests that t he /^-dependence of whistler heat flux 
instability (|Garv et al.l Il994h is inconsistent with 
the data, however a magnetosonic instability may 
be consistent in the collisionless regime. We also 
infer a /3 e -dependance to the temperature profile 
index a, which may indicate an additional en- 
ergy transport process, or be a proxy for another 
plasma parameter (e.g. Mach number, collisional- 
ity, etc.). 

The observation that the electron heat flux re- 
mains 'classical' to Knudsen numbers as large as 
Kt ~ 0.3 may have implications for solar wind 
models which transition from fluid to exospheric 
domains ( viz. Echim et al.l2011 and references therein! ) 
and for overall heat transport in the corona. If 
the electron temperature gradient is only a weak 
function of radial distance (in the free solar wind), 
then Xfp will be smaller in the inner heliosphere, 
and it may be that most of the solar wind lies in 
the SH regime there. 

Electron thermal conduction, in both the col- 
lisional Spitzer-Harm and collisionless limits, is 
energetically and dynamically important in other 
low-collisionality astrophysical plasmas, including, 
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e.g., the hot intracluster medium in gal axy clus- 
ters ( e.g. Bertschinger and Meiksinlll986h . the hot 
interstellar medium in galaxies (Cowie & McKee 
1977), and some accret ion disks around neutro n 
stars and black holes ( e.g. Sharma et all 12006). 
These astrophysical plasmas are probably charac- 
terized by j3 e ~ 1 — 10, not too dissimilar from a 
large fraction of the epochs of in situ solar wind 
data used here (Table 1). Our results suggest that 
parallel thermal conduction is likely to be compa- 
rable to the Spitzer-Harm and/or saturated values 
in these systems. In more detail, the results in Fig. 
3 and Table 1 could be used for modeling paral- 
lel thermal conduction in other astrophysical con- 
texts. Of course, the global magnetic field geom- 
etry determines in part how the parallel thermal 
conductivity translates into large-scale redistribu- 
tion of heat (e.g., Chandran & Maron 2004), and 
the field geometry is significantly more uncertain 
in these other astrophysical contexts. 

This work was supported in part by NSF grant 
AGS 0962726 (SHINE) and NASA HTP grant 
NNX11AJ37G. 
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